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Abstract 

We develop a formalism to study non-Gaussianity in both curvature and isocur- 
vature perturbations. It is shown that non-Gaussianity in the isocurvature pertur- 
bation between dark matter and photons leaves distinct signatures in the CMB tem- 
perature fluctuations, which may be confirmed in future experiments, or possibly, 
even in the currently available observational data. As an explicit example, we con- 
sider the QCD axion and show that it can actually induce sizable non-Gaussianity 
for the infiationary scale, i^inf = 0(10^ - 10^^) GeV. 



1 Introduction 



The accumulating observational data, especially the WMAP observation of the cosmic mi- 
crowave background (CMB) [Ij, provided significant support for the infiationary paradigm. 
The results of these measurements are consistent with nearly scale-invariant, adiabatic 
and Gaussian primordial density perturbations, known as the standard lore in the simple 
class of infiation models. 

A possible detection of the deviation from the above properties will enable us to further 
constrain infiation models. The scalar spectral index of the power spectrum is constrained 
as Tis = 0.963lo!i5 68% C.L. [Ij, which already excludes some infiation models. No 
significant isocurvature component has been detected so far, and the current constraint 
on the ratio of the amplitudes of the isocurvature and curvature perturbations reads, 
l^/Cl ^ 0.3 m Ej. Recently, Yadav and Wandelt claimed an evidence of the significant 
non-Gaussianity in the CMB anisotropy data. Using the non-linearlity parameter /nl to 
be defined in the next section, their result is written as 47 < /nl < 127 at 95% C.L. [3j. 
On the other hand, the latest WMAP five-year result is consistent with the vanishing 
non-Gaussianity: —9 < /nl < HI at 95% C.L., including /nl = 0. Interestingly, 
however, the likelihood distribution of the WMAP result is biased toward positive values 
of /nl- Also there are some other studies searching for the non-Gaussianity [4J, and it is 
not settled yet whether the non-Gaussianity exists. At the present stage, therefore, it is 
fair to say that the observations are consistent with the nearly scale invariant and pure 
adiabatic perturbations with Gaussian statistics, while there is a hint of non-Gaussianity 
at the two sigma level. 

The standard lore on infiation is based on a simple but crude assumption that it is 
only the infiaton that acquires sizable quantum fiuctuations during infiation. Its apparent 
success, however, does not necessarily mean that such a non-trivial condition is commonly 
met in the landscape of the infiation theory. In fact, there are many fiat directions in a 
supersymmetric (SUSY) theory and the string theory. If some of them are light during 
infiation, they acquire quantum fiuctuations, which may result in slight deviation from 
the standard lore. 

One promising candidate is provided by the theory with a Peccei-Quinn (PQ) symme- 
try, which is introduced in order to solve the strong CP problem in the quantum chromody- 
namics (QCD) [5[|6jJ. There appears a pseudo-Nambu-Goldstone boson called axion associ- 
ated with the spontaneous breakdown of the PQ symmetry. The axion is a light scalar field 
and contributes to the cold dark matter (CDM) of the universe [7J. In particular, the axion 
can have a large isocurvature perturbation 0,[9]. As for non-Gaussianity, it is known that 
the slow-roll infiation generally predicts a negligible amount of non-Gaussianity, /nl = 
0(e, 77) [TOl [m, [12l [13]. Here e and rj are the slow-roll parameters, which must be smaller 
than unity for the slow-roll infiation to last long enough. In the curvaton [HI [151 EB] 
and/or ungaussiton [T7] scenarios, however, there are light scalars in addition to the in- 
fiaton, which can generate sizable non-Gaussianity [IHlIiniEOliaiESlESllHIIIl- As we 

^ Here we have quoted the value of /nl^^ since we are interested in non-Gaussianity of the local type 
in this paper. 
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will see, the axion can also induce sizable non-Gaussianity. 

In this paper we point out that, if an isocurvature component possesses some amount 
of non-Gaussianity, it is transferred to the non-Gaussianity of the curvature perturbation, 
resulting in a possibly large value of /nl- If there are no other light scalar fields than the 
infiaton, the resultant density perturbations are necessarily adiabatic and almost Gaus- 
sian. As mentioned before, this may not be the case in the presence of many fiat directions. 
Suppose that there is a light scalar that acquires quantum fiuctuations during infiation. 
Then its fiuctuations produce isocurvature component. If the scalar decays into radiation, 
the isocurvature perturbation is converted into the adiabatic one. This is exactly what 
occurs in the curvaton and/or ungaussiton scenarios. However, as far as non-Gaussianity 
is concerned, the light scalar having large fiuctuations needs not decay. Even if such 
a light field is not responsible for the total curvature perturbation, it can still provide a 
source of large non-Gaussianity. This interesting possibility was noted in Refs. p!5ll25ll26] . 
In this paper we have systematically studied the non-Gaussainity from the isocurvature 
perturbations and how it exhibits itself in the CMB anisotropy. Interestingly enough, we 
have found that the resultant non-Gaussianity induced by the isocurvature perturbations 
has distinctive signatures in the CMB, which should be distinguished from that in the 
curvaton and ungaussiton scenarios. 

This paper is organized as follows. In Sec. O a general formalism to study non- 
Gaussianity including isocurvature perturbations is presented. In Sec. [3l we compute the 
bispectrum of the temperature fiuctuations arising from non-Gaussianity in the isocurva- 
ture perturbations. In Sec. H] the formalism is applied to the case of the axion and it is 
shown that the axion can induce large non-Gaussianity while leaving a certain amount of 
the CDM isocurvature perturbation. Sec. [5] is devoted to discussion and conclusions. 

2 Non-linear isocurvature perturbation 
2.1 Definition of the isocurvature perturbation 

Let us consider cosmological perturbations of multicomponent fluids labeled by i = 
1, . . . , n. We assume that the density perturbations originate from fluctuations of scalar 
flelds generated during inflation. 
We write the spacetime metric as 



where Af is the lapse function, the shift vector, the spatial metric, a{t) the back- 
ground scale factor, and i/j the curvature perturbation. On sufliciently large spatial scales, 
the curvature perturbation ip on an arbitrary slicing at t = t/ is expressed by 



ds^ = -N^dt^ + a^{t)e^^-i,j {dx' + I3'dt) {dx' + [3'dt) 



(1) 



ipitf.x) = N{tf,ti : x) - log 



(2) 



a{ti) 
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where the initial shcing at t = is chosen in such a way that the curvature perturbations 
vanish (flat shcing). Here N{tf^ti : x) is the local e-folding number, given by the integral 
of the local expansion along the worldline x = const, from t = ti to t = tf. 

We denote by ( the curvature perturbation ijj evaluated on the slice where the total 
energy density is spatially uniform (uniform-density slicing). In a similar fashion, we also 
introduce Q to denote the curvature perturbation on the slice where pi is uniform {6pi = 
slicing). Then, from Eq. ([2j), ("^ is related to ( by the gauge transformation 

Ci = C + ANi, (3) 

where AA^^ is the e-folding number measured from the uniform-density slicing to the Spi = 
slicing, both slicings corresponding to the same background time. If each component of 
the fluids does not exchange its energy with the others, Q are known to remain constant 
for the scales larger than the horizon [27j. 

Let us deflne 6pi as the perturbation of pi on the uniform-density slicing. That is, 

p^{N,x) = p,Q + Sp,{N,x), (4) 

where pio is the energy density of the i-th fluid in the background spacetime, and = 
N{x) deflnes the uniform density slicing. Then AA^^ is related to Spi by the following 
equation, 

p,{N + AA^„ x) = p,{N, x) - Sp,{N, f), (5) 

where the l.h.s and r.h.s are evaluated on the 6pi = slicing and on the uniform-density 
slicing, respectively. Assuming Spi/pio <C 1 and AA^^ <C 1, we can solve this equation 
with respect to AA^^ up to the second order in Spf. 

PiO ^PiO \ PiO J PiO 

where the prime denotes the derivative with respect to Hence Q can be written as 

Spi PiO fSPiV , Spi^P'i 



PiO ^PiO \ PiO J PiO 

We define the (non-Unear) isocurvature perturbation between the i-th fluid and the 
j-th one as [28] 

5,, =3(0-0)- (8) 
Using Eq. ([71), Sij can be written as 



Sij — 3 



P'io^ P'jo 2p^oVPio/ ^ '^p'joKP'joJ ^ p'io P'jo 



(9) 



^ If the z-th fluid has vanishing homogenous value, i.e., if it is produced predominantly by the quantum 
fluctuations, 5pi/ piQ as well as AA^^ is no longer small. In the example of axion which we discuss later, 
this problem can be avoided by considering the density contrast of the total CDM sector. 
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If we neglect the second order terms, 5^^^ reduces to the weU-known form. If the i-th 
fluid fluctuates in the same way as the j-th one, i.e., AA^^ = ANj = 0, the isocurvature 
perturbation between the two, Sij^ vanishes. AU the isocurvature perturbations vanish if 
there is only the adiabatic perturbation, that is, if all AA^^ vanish. 

We assume that the density perturbations originate from the fluctuations of light scalar 
fields during inflation. Then Spi can be expanded as[f| 



Pi,aC 



+ 



(10) 



where is the quantum fluctuation of a light scalar 0^ on the initial flat slicing at 
t = ti. We choose the initial time ti slightly after the cosmological scales of interest exit 
the Hubble horizon, since the above formulation is valid for the superhorizon modes. We 
assume that the scalar flelds, {(/>a}, have quadratic potential and behave like free flelds 
and do not have any sizable interactions during inflation. In particular, their masses are 
assumed to be lighter than the Hubble parameter during inflation. In this case, the higher 
order terms in Eq. (fTOl) are safely neglected. Then, to a good approximation, 5(/)a is given 
by the Gaussian variable. In general, the coeflicients that appear in the right-hand side 
of Eq. (fTOl) depend on the slicing on which they are evaluated. When evaluating the 
coefl&cients, we need to choose an appropriate uniform-density slicing. For example, if 
Pi denotes the energy density of the axion, those coefl&cients are easily evaluated on the 
uniform density slicing when the axion starts to oscillate. If pi is the energy density of a 
particle produced by the decay of a scalar fleld, those coeflicients include the information 
from the onset of the flled oscillation to its decay. Thus case-by-case calculations are 
required. Substituting Eq. (fTOl) into Eq. ([9]), Sij can be written in the form 



Sij 



6" + -Sij^ab^ 



(11) 



with 



2 I Pi,a _ Pj,a 



= -3 



+6 



PiO PjO 

Pi,ab Pj,cib 
PiO PjO 
Pi^aPi^b Pj,aPj^b 



p% 



PiOpi,aPi,b PjOPj,aPj,b 



P iO 



PjO 



(12) 



(13) 



For simplicity, we assume that the masses of are negligible, and the fluctuations 
are independent to each other. Then the correlation functions are given by the following 
form, 

{ScPlScPl) = {27:fP,^{h)8{h + h)5^' (14) 



Note that {5pi] are subject to a constraint 5 pi = on the uniform density shcing since. 
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with 

PsM ^ (15) 

where k denotes the comoving wavenumber, and Hi^f is the Hubble parameter during 
inflation. For later use, we also deflne the following: 



2.2 Bispectrum of the isocurvature perturbations 

We deflne the power spectrum and bispectrum of Sij as 

(S.fkAfJ = {2nfPs^,{h)Sih + k2), (17) 

and 

(S^fkAjkAjh) ^ i2^fBs^J{ku k2, k,)S{h + k2 + h). (18) 

Here and in what follows no summation is taken over the indices i and while we sum 
over the repeated indices a^b^c^.... Using ( |TTI ) and ( |T4l ). the power spectrum can be 
expressed as 

1 f dP'k' ^ ^ 

Psij{k) = Sij^aSij^aPs<p{k) + -Sij^abSij^ab J j^^Ps<p{k') Ps<p{\k - \) , (19) 

where we have regarded Scj)^ as a Gaussian variable. After performing the integration, we 
obtain 

Psij{k) = [Sij^aSij^a + Sij^abSij^ab^^Sct>^n{kL)] Pscf,{k), (20) 

where we have introduced an infrared cutoflF L that is taken to be of order of the present 
Hubble horizon scale [29l [26[ |30] . Similarly, the bispectrum can be written as 

Bsij{ki, ^2, ks) = Sij^aSij^bSij^ab [P6ct>{ki)P6cf>{k2) + (2 perms)] 

+ / ||p,,(.')P..(|J: - -t'M + (21) 

In the squeezed conflguration in which one of the three wavenumbers is much smaller than 
the other two (e.g. fci <C /c2, ^3), it is approximately given by 

Bsij{ki^ /C2, ^3) — \Sij^aSij^bSij^ab + Sij^abSij^bcSij^ca^Sct) ^^(^6^)] 

X [Ps^{ki)Ps^{k2) + P54>{k2)Ps^ik3) + P54,{k3)Ps^{ki)] , (22) 

where kb = min{/ci, ^2,^3}- 

Let us define the non-hearity parameter of the isocurvature perturbations, fsij, as 

Bsijiki, k2, ks) = fsiAPsij{ki)Psij{k2) + Psij{k2)Psij{kz) + Psij{ki)Psij{ki)]- (23) 
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We can see that fsij is not very sensitive to the wavenumbers. If Sij is dominated by the 
hnear terms in (see (fTTI)). fsij becomes independent of the wavenumbers, and given 
by 

r ^ Sij^JjSij^cSij^Jjc fOA\ 

jSij — —rZ - 72" 5 1^4 j 

[^ij,a^ij,a) 

for generic configurations of the wavenumbers. Even if the quadratic part dominates, i.e., 
Sij 1/2 Sij^ah5(j)^5(j)^ ^ its dependence is only logarithmic in the squeezed configuration: 

1 Sij^abSij^bcSij^ca /OP;\ 

where we have approximated as ln(/c^L) lii{ki,L) for i = 1,2,3. For a generic config- 
uration, the dependence may become more involved. Nevertheless, we expect that such 
dependence is also mild for the scales of interest, based on the dimensional arguments. 



2.3 "/nl" from fs^j 

In many literatures, the non-linearity parameter /nl is used to measure the non-Gaussianity 
of the adiabatic perturbations (for example, see Ref. [31j and references therein). We 
adopt the following conventional definition of /nl, 

iCkSk.CkJ = lui^^fSih + h + h) [PdkdPdh) + 2 perms.] . (26) 

Now we would like to relate /nl to fsij defined by Eq. ( l23l ). This is a non-trivial task 
since we are considering the non-Gaussianity of the isocurvature perturbations. To be 
definite, we hereafter consider the CDM isocurvature perturbations. It can be extended 
to the other types of the isocurvature perturbations in a similar way. 

We write the adiabatic perturbations originating from the infiaton fiuctuations as C^^. 
Using the formula (Ej), it can be expressed as 

r' = N^H, (27) 

where 5(j) denotes the fiuctuation of the infiaton 0, and A^^^ is the derivative of the local e- 
folding number with respect to (j). We assume that primordial infiation does not generate 
large non-Gaussianity. In fact, it was shown that the non-linearlity parameter generated 
during slow-roll infiation is at most of the order of the slow-roll parameters, and hence 
below the sensitivity of the Planck satellite | lU [121 [13] . Then the three-point function of 
Q^^ is approximately given by 

Let us now define the CDM isocurvature perturbation in the universe after the reheat- 
ing as 

S = 3(CCDM-Cr-), (29) 
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where CcDM(r) is the curvature perturbation on the shcing where the energy density of 
the CDM(radiation) is spatiaUy uniform. We assume that the CDM is always decoupled 
from the radiation, so that Ccdm as well as (r are time-independent. Note that ( is not 
necessarily conserved in the presence of the isocurvature perturbation. 

When the universe is dominated by the radiation, the curvature perturbation on the 
super-horizon scales is given by ( = (r- We assume that the curvature perturbation at 
that time is originated solely from the inflaton, i.e., (r = C^^^- the matter dominated 
era, we have ( = Ccdm- Hence ( in the matter dominated era can be written as 

1 

C = -\ — (matter dominated era), (30) 
3 

where we have assumed that the curvature perturbation mainly comes from the inflaton 
and the other flelds contribute only to the isocurvature perturbations. Note that this 
relation holds to any orders in the perturbative expansion. 

We assume that the isocurvature perturbation is uncorrelated with the primordial 
curvature perturbation, i.e., 

(Qf^fe) = 0- (31) 

The three-point function of the curvature perturbation is then evaluated as 

i^kSk2^k3^ ^ '^^^ki^k2^k3y (22) 

From (l23l) . (I26l) and (l32l) . /nl is related to fs as follows, 

6 ^ 1 Ps{ki)Ps{k2) + (2 perms) 

g/NL 27^^ Pdh)P^{k,) + (2 perms) ^^^^ 

= ^^fs [S^aS^a + S^abS^ab^^ ln(/cL)] ^ . 

It should be noted that the above relation between /nl and fs is valid only for the large 
scales which enter the horizon after the matter-radiation equality. This however helps us 
to get a feeling of the non-Gaussianity produced from the isocurvature perturbation. 



3 CMB Temperature Fluctuations 

In this section, we calculate how the non-Gaussianity of the isocurvature perturbation 
exhibits itself in the CMB anisotropy, following the notations used in Ref . [HH [32j . In 
particular, it contributes to the bispectrum of the CMB temperature fluctuations, which 
may be observable in the future observations. 

We introduce spherical harmonic coeflicients of the temperature anisotropy arising 
from the CDM isocurvature perturbations as 

<^im = J dri Y^M. (34) 
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In order to relate the primordial fluctuations to the CMB temperature anisotropy, one 
needs to multiply the transfer function. Precisely speaking, one has to use a non-linear 
version of the transfer function, which also induces a certain amount of non-Gaussianity. 
Based on the dimensional grounds, however, such secondary non-Gaussianity is expected 
to be much smaller than the value currently hinted by the observation. Since we are 
interested in the relatively large primordial non-Gaussian features in the CMB anisotropy, 
we can neglect the intrinsic non-linear property in the transfer function. We therefore use 
the linear transfer function, g^^\k)^ deflned by 

ef°)(fc) = g^ft^m^ (35) 
where B^^^^^(fc) is the multipole moment of CMB temperature anisotropy: 

^ = /(§i:''p^-i)er'm(£^«). (36) 

Here P^'s are the Legendre polynomials. Using Eqs. f l35l ) and f l36l ). a^^^^ can be written 

as 

(27r)= 



j ^,9'T^\k)YL{i)Sn- (37) 



The anglar power spectrum of a^^^^ is deflned by 

/^(iso) (iS0)*\ _ ^(iso)r r /00\ 

Using (l37l) . we obtain 

Cf^^ = - / dk k\g^f,^\k)) Ps{k). (39) 
^ Jo 

Here and in what follows, we use A; = to show the lower limit of the integration interval, 
although it is set to be the infrared cutoflF in the actual calculation. The angular 
bispectrum of a^^^"^ is deflned by 

/^(iso) (iso) (iso) \ _ p(iso)mim2m3 fAr^\ 

Statistical isotropy divides the angular bispectrum into the following form, 

o(iso)mim2m3 _ ^mim2m3 l^so) / . \ 

Here G^l^^^^ = J dn Y^^^^(n)Y^2m2(^)^3m3(^) (Gaunt integral) and ^^^^^2^3 reduced 
bispectrum, on which we will focus in the following. 
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Substitutiing (1371) into (HOl) . we obtain 

o roo roo roo roo 

teL = -3 drr' dk, / dk^ kl / dk, kj 
^ Jo Jo Jo Jo 

x^T? (^Oi^i (^i^)^r? (^2)i^2 {k2r)gT7^ {h)3h ik^r)Bs{ki, fca, fcs) 
= y (irr^y dki klg!^l"^{ki)je^{kir)Ps{ki) j dk2 k^g'^^'^J {k2)je^{k2r)Ps{k2) 

POO 

X / rffca fc^^Jf (fc3)j,3(fc3r) /^(fci, A;2, ^3) + (2 perms), (42) 
Jo 

where ji{x) is the spherical Bessel function, and we have used Eq. ( I23l) in the last equal- 
ity. From the discussion below Eq. (I22l) . we have seen that fs depends on the three 
wavenumbers at most logarithmically, i.e. the dependence is rather weak. Therefore, we 
approximate fs as the constant and write the bispectrum as 

tSe. - fs I dr (6^:^)^ (0^1(0 + 2 perms.) , (43) 

where 6^^^^ and b^^-^j are defined by 

blT\r) = - / dkk'g!f;\k)Mkr)Ps{k), (44) 

o poo 

r(iso)/ X _ ^ / 7 7, 7 2^(iS0) 



TT 



= - ciA:A:^^;^r(A:)j,(A:r). (45) 



In a similar way, we can define and W^i^i^ as the angular power spectrum and the 
reduced bispectrum of the total temperature anisotropy including both the adiabatic and 
the isocurvature contributions. We define the non-linearity parameter as 



6 

^hhh = -^Inl J dr {bi^i,{r)bu2{r)b^Usir) + 2 perms.) 



(46) 



with 



bu{r)^- dkk'g^^t\k)u{kT)P^{k), (47) 
Jo 

r) poo 

b^Ur)^- dkk^g^^f{k)j^{kT), (48) 
Jo 

where gTt\k) is the transfer function for the adiabatic perturbations. Note that it is 
f^^ that is directly related to the CMB observations If only the adiabatic perturbation 
exists, /^j^ coincides with /nl defined by Eq. f l26l ). However, the relation between f^^ 



^ A rigorous procedure to constrain non-Gaussianity from isocurvature perturbations using observa- 
tional data can be done in a way demonstrated in Ref. [33j . Here we simply use as a representative 
value which characterizes the non-Gaussianity in the CMB anisotropy. 
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Non-linearity parameter 


Related to 


Definition 


fs 
/nl 


3-point function of isocurvature perturbation 
3-point function of curvature perturbation 
3-point function of temperature perturbation 


Eq.([23l) 
Eq.(|46]) 



Table 1: Non-linearity parameters. 



and /nl gets rather involved when the isocurvature perturbation mainly contributes to 
the bispectrum while the power spectrum is dominated by the primordial adiabatic con- 
tribution, i.e., Ci ^ (7^^^^^^ and W^i^i^ ^hhh' particular, it should be noted that /^^^ 
sensitively depends on (^1,^2,^3)- 

Table d] summarizes the non-linearity parameters which we have defined so far, fs^ /nl 
and /nl^. Given a model, fs is easily calculated by Eqs. (IMI) and (l25l) . Once we know 
fs we can obtain /nl through the relation (l33l) . But the most relevant quantity directly 
related to the CMB observations is f^^ and we evaluate it numerically in the following. 



3.1 Sachs- Wolfe approximation 

For the low multipoles, typically smaller than ^ ~ 10, the temperature anisotropy comes 
mainly from the Sachs- Wolfe effect. 

From this equation, the transfer function in the Sachs- Wolfe regime can be written as 

9?;\k) = -jjeikn), (50) 

where is the comoving distance to the last scattering surface from us. Then the reduced 
bispectrum becomes 

,(iso) 5 /^(iso)^(iso) 



& - -T^fs Crcr + 2 perms. . (51) 



As we will see in the next section, the above expression f lSTl ) is not that precise since 6nl^(^) 
has non-negligible contributions from smaller scales beyond the Sachs- Wolfe plateau. It is 
still useful however to understand the large amplitude of the bispectrum from isocurvature 
perturbations with non-Gaussianity. Here we continue with this approximation and derive 
a relation between f^^ and /nl which is valid up to a 0(1) numerical factor, and leave 
detailed discussions for the next subsection. 

Under this approximation, f^^ defined in Eq. ( H6l) is expressed as 

hfM = -6/nl + 2 perms.) . (52) 

We require that the adiabatic perturbations dominate the power spectrum, i.e., 
^^adi)^ When only the adiabatic perturbations exist, f^^ exactly coincides with /nl [31j . 
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However, in the presence of the isocurvature perturbations, f^^ is different from /nl- If 
the bispectrum is 
can be written as 



the bispectrum is dominated by the isocurvature perturbation, i.e., b^^^^^^ ^hhh^ -^nl^ 



5 C^r'c1r'+ 2 pern. 
-12cr'cr' + 2perms/"- ' ' 

Combining this relation with Eq. (I33l) yields 

fil ~ 216/nl. (54) 

This is a quite important result. When the non-Gaussianity comes from the isocurvature 
perturbations, the non-Gaussian features appearing in the CMB anisotropy is greatly 
enhanced in the Sachs- Wolfe plateau. This is because the isocurvature perturbations 
make more contributions to the CMB power spectrum at low multipoles than the adiabatic 
perturbations. As we will see below, full treatment of the transfer functions beyond the 
Sachs- Wolfe plateau will change the numerical factor 216 in Eq. (I5^ to about 100 . 



3.2 Acoustic scales 

In this section we present more detailed discussion on the bispectrum from the non- 
Gaussian isocurvature perturbations, especially focusing on their differences from the 
adiabatic perturbations. To study the features in the bispectrum at small angular scales 
beyond the Sachs- Wolfe plateau, we have numerically calculated 6L^(r), 6nl^(^) and the re- 
duced bispectrum b£^£^£^ using transfer functions gTi{k) from the CAMB code [3lj. Through- 
out this section we adopt the flat SCDM model and assume a set of cosmological pa- 
rameters {Qb = 0.05, Qc = 0.95, h = 0.5), where ^b(c) is the density parameter of the 
baryon(CDM), and h is the Hubble parameter in units of lOOkm/sec/Mpc. For simplic- 
ity, we neglect the tilt of the power spectra, Ps{k) and P({k) in Eq. (I33l) . Then both 
Ps{k) and P({k) are proportional to k~^. 

In Fig. [1] we show the numerical results of 6L^(r) and 6nl^(^)- If we take r largely 
diflFerent from = tq — r*, which is the comoving distance from us at r = tq to the 
last scattering surface at r = r*, both and 6nl^ get suppressed. This is because 
dTiik) ^ ji{kr^) and ji{kr) in the integrant of Eqs. (I44II45I) would oscillate with diflFerent 
frequencies, making contributions in a destructive way. Therefore the signature of the 
primordial non-Gaussianity mostly comes from r r^. In Fig. [1] we have taken several 
values of r around r^. 

Let us flrst consider bu{r) shown in the upper panels of Fig. [H We notice that fc^T^ 
is roughly twice as large as b^^^'' in the amplitude at large angular scales {£ < 10) for 
any values of r r^^. This behavior is similar to the angular power spectra C^. It can be 
easily understood by noting that the Sachs- Wolfe effect leads to 

g!f;\k)c.2gi^^Hk). (55) 



11 




Figure 1: 6L^(r) (top) and 6nl^(^) (bottom) from numerical calculation. We show the 
cases of isocurvature initial conditions (solid red line) and adiabatic conditions (dashed 
green line). Here we have assumed Ps{k) = Pc{k). From left to right panels, r is set to 
be r = To — 0.5r*, tq — and tq — 1.5r*, respectively. Note that is a dimensionless 
quantity, while 6nl^ has the dimensionality of (length)"^. 



At smaller angular scales 6L^(r)'s represent the acoustic oscillations. Note that the phase 
of oscillations are different by 7r/2 between the isocurvature and adiabatic perturbations. 
This is another similarity of 6L^(r) with C^. 

On the other hand, the situation with is slightly different. Athough 6nl^'s also give 
similar flat spectra at the large angular scales, they have some differences from 6L^(r). 
One of them is that the ratio of the amplitudes of 6nl^ for isocurvature and adiabatic 
initial conditions differs with r. This is because 6nl^(^) receives more contribution from 
smaller scales due to the absence of in the denominator of the integrant in Eq (Il5l) . 
compared with 6L^(r). Since the perturbations in photon fluid becomes smaller at large 
k for isocurvature perturbations, the b^^-^^j is not as large as 26^^^^ This changes Eq. ( I54l) 
obtained by using the approximated transfer function in the Sachs- Wolfe regime Eq. (I50l) . 
We will discuss this issue below. 
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Figure 2: The reduced bispectra ^2(^2 + 1)4(4 + (^^p) and + 

1)4(4 + l)6^;^J^^3^^/(27r)2 (bottom). To avoid complexity, we have fixed (^i,4) = (9, H) 
(left), (99, 101) (middle), (199,201) (right) and varied 4. The solid red line and dashed 
green line correspond to the cases with isocurvature and adiabatic initial conditions, 
separately. The unobservable multipoles are shown as shaded region. We have set fs = l. 



Among the three terms in Eq. f l43l ). we denote the two of them by 



7L,L,NL 



= J drr^bu,{r)bu2{r)bNLh{r)- 



(56) 



In Fig. O we show the bispectrum b^^/J^^^ and b^^j'^^^ for the isocurvature and adiabatic 
perturbations. One can see that the amplitude of the bispectrum of the isocurvature 
perturbations is enhanced at large angular scales (^i,4,4 ^ 10), compared with the 
adiabatic perturbations. Our numerical calculations give the ratio at large scales as 



7 (adi) 



4. 



(57) 



If the transfer function g^^\k) in Eq. f ISOl ) were valid at smaller scales beyond the Sachs- 
Wolfe plateau, the right hand side of Eq f l57l ) would be 8, as expected from Eqs. f l49l ) and 
f l55l ). However as have mentioned above, 6nl^(^) receives contributions from smalle scales 
{k > I/a/St*) in a destructive way, and this roughly halves the the amplitude of the bis- 
pectrum from isocurvature perturbations compared to that from adiabatic perturbations. 
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Figure 3: /nl^/Znl is plotted as a function of ^3 with various sets of (^1,^2)- Only 
observable multipoles — ^2! ^ ^3 ^ ^1 + ^2) are shown. 



At acoustic region (^1,^2,^3 — 200), the amplitude of ^^^72^3 becomes much suppressed 
since CMB temperature anisotropics are small at acoustic region for isocurvature initial 
perturbations. 

In Fig. [3] we plot /nl^//nl as a function of is for various sets of (^1,^2)- We can see 
that our estimate of at large angular scales in Eq. (I5^ is corrected and approximately 
given by 

f^l ^ 100/nl. (58) 

However the amplitude of the bispectrum at large angular scales is still enhanced for 
non-Gaussian isocurvature perturbations compared with adiabatic perturbations. In the 
adiabatic case, the largest signals of the primordial non-Gaussianity come from the acous- 
tic regions. On the other hand, in the isocurvature case, the signals are concentrated in 
the large angular scales. This remarkable difference in -^-dependence of ^^^^72^3 ^hhh 
will help us distinguish the isocurvature non-Gaussianity from the adiabatic one. 

So far, we have assumed SCDM universe. Here we give some comments on how the 
above results will change in the standard ACDM universe. Assuming the flat universe, 
the amount of CDM is smaller {Qc — 0.21 pLj) in the ACDM universe. For the adiabatic 
initial perturbations, this does not cause much diflFerence on the CMB anisotropics at large 
angular scales, except for the small late-time Integrated Sachs- Wolfe effect from nonzero 
f^A. On the other hand, for the isocurvature initial perturbations, the anisotropics at 
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large angular scales will be suppressed. This is because the universe is not completely 
matter-dominated at recombination and therefore the curvature perturbations on large 
scales are not entirely generated from isocurvature perturbations in CDM (Eq. (I30l) is 
not a very good approximation in the ACDM model). At large angular scales, the CMB 
anisotropics with the isocurvature initial perturbations are, however, still larger than for 
adiabatic ones and our discussion above is basically valid even in the ACDM universe. 
When constraining the isocurvature non-Gaussianity by using the future CMB data, we 
should take more realistic ACDM model, but there will be no fundamental difference. 



4 Application to the axion 

In this section we apply our formulation to the axion as a concrete example. The axion, 
a, is a pseudo Nambu-Goldstone boson associated with the spontaneous breaking of the 
PQ symmetry. Let us denote the breaking scale by F^, whose magnitude is constrained 
from various experiments, astrophysical and cosmological considerations. The most strict 
lower bound on Fa comes from the observation that the duration of the neutrino burst 
in SN1987A lasted for '-^lO seconds. In order to prevent too fast cooling by the axion 
emission. Fa ^ 10^° GeV is required [35j. On the other hand, the upper bound is provided 
by the cosmological argument. The axion obtains a tiny mass after the QCD phase 
transition due to the anomaly effect which explicitly breaks the PQ symmetry. The axion 
begins to oscillate coherently after that. Since the lifetime of the axion is very long, it 
survives until now and contributes to DM of the universe. The abundance is estimated 
as [36] 



1.7 \ 1-18 

1012 G^ 



0-2(t7§^) for Fae>H,^,/27i, 



Fa \-'-'' ( H,^,/2^ Y 



(59) 



where 9 denotes the initial misalignment angle of the axion. Imposing Qah < 
0.11 [Ij, we obtain an upper bound on the PQ scale as Fa ^ ^"^ '^lO^^ GeV. Thus the ratio 
of the axion abundance to the total dark matter abundance r = f^a/^CDM is calculated 
as 

-2f Fa V^'^\ ^* Vf 0.11 



where we have defined 



a* = max <{ F«(9, } . (61) 



If the PQ symmetry is already broken before or during infiation and if it is never 
restored after infiation, the axion has unsuppressed quantum fiuctuations during infiation 
because it remains practically massless during infiation. Since the axion contributes to 
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some fraction of DM, such an axionic isocurvature fluctuation is converted to the CDM 
isocurvature fluctuation. Thus the axion is a plausible candidate for generating the non- 
Gaussianity from the isocurvature perturbation. 

Now let us estimate the magnitude of the axionic isocurvature fluctuation and the 
resultant non-Gaussianity. We assume that the inflaton itself does not generate non- 
Gaussianity, and that only the axion has an isocurvature fluctuation. The axion acquires 
a quantum fluctuation during inflation given by 

ALW ^ (62) 

The observationally relevant quantity is the CDM isocurvature perturbation, rather than 
the axionic isocurvature perturbation itself. Using Eq. ([9]), the CDM isocurvature per- 
turbation S is given by 

Sc^ ^ (63) 

PmO + PaO 

where Pmo denotes the dark matter abundance other than the axion. In this equation, 
5 Pa is evaluated on the uniform density slicing. Since S is conserved quantity as long as 
the scales of interest are sufliciently large, we can evaluate it when the axion starts to 
oscillate. When the axion starts to oscillate, the universe is dominated by the radiation, 
and therefore we can safely neglect the density fluctuation of the radiation, 5pr/ Pr- 

As will become clear later, if one imposes the current constraint on the isocurvature 
perturbation, large non-Gaussianity is generated only for r <C 1. Using r <C 1, we can 
approximate S as 



S 



2- 



2ai6a / 6a\ 



(64) 



where = Fa9 denotes the classical deviation from the potential minimum. The flrst 
term in f iGTl ) dominates over the second one when the classical deviation from the potential 
minimum overcomes the amplitude of the quantum fluctuation. In the opposite case, i.e., 
the isocurvature fluctuation is dominated by the second term in (IMl) . the whole dynamics 
of the axion is controlled by the quantum fluctuation generated during inflation. This is 
similar to the "ungaussiton" scenario [17j: the axion is predominantly produced by the 
quantum fluctuations, giving the non-Gaussianity to the density fluctuations, while its 
contribution to the total curvature perturbation is negligibly small. 

For later convenience, the power spectrum of the isocurvature perturbation Ps{k) is 
deflned through Eq. (fTTI) . 

{Sj:^Sj:^) = {2nmh + h)Ps{h). (65) 

Also note that the isocurvature perturbation is uncorrelated with the primordial curvature 
perturbation in the case of the axion. 
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It is straitforward to calculate the non-linearity parameter /nl defined by Eq. (I26l) . 
First note that from Eqs. IHM and (l25l) fs is calculated as 



fs 



1 

2r 



for Fae>Hi,i/27i 



1/2 



2r \AsaJ HkL) V^s(fc); 

(67) 

where 2.5 x 10 ^ is the WMAP normalization of the curvature perturbation. Then 

Eq. f l33l ) tells us that the non-linearity parameter is given by 

,,, = ^f§WV. (68) 



324r V^c(^) 

if the classical deviation overcomes the quantum fiuctuation [FaO > i7inf/27r). Since 
the WMAP five-year results give a contraint Ps/Pq % 0.1, a small value of r < 10~^ is 
necessary for large non-Gaussianity, /nl ^ But there is an upper bound on the level 

of non-Gaussianity coming from CDM isocurvature perturbation. In order to maximize 
the value of /nl, the isocurvature perturbation must also be large. However, since S is 
limited as < r, /nl is maximized at 2r ^ 0.3A^ in order to saturate the isocurvature 
bound and this gives a strict upper bound as /nl ^ 20. This is explicitly shown in the 
case 5a ^ Fa9^ where the quantum fiuctuation dominates the axion dynamics, giving 
S r. In this case the non-linearlity parameter is estimated as 



^S/SW)''^l„(fcL)]-v^ 




(69) 



Note that in this regime the parameter dependence of /nl is all compressed in the informa- 
tion of the magnitude of the isocurvature perturbation. Thus non-Gaussianity parameter 
/nl is solely bounded by the isocurvature constraint in this regime. In other words, /nl 
is maximized when the isocurvature contribution saturates the allowed maximum value 
and the bound /nl ^ 20 does not depend on other model parameters. It may also be 
useful to give a full expression for /nl, 



(70) 



whose limiting behavior approaches to the expressions given above. 

Above results can be understood in a simple way. The following rough estimations 
may be useful because they give a correct parameter dependence. When 6a <C Fa9 the 
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Hint [GeV] 



Figure 4: Contours of /nl = 1, 10 and 100 for Fa = 10^° GeV. Gray shaded region 
is excluded from isocurvature constraint. In the blue region the PQ symmetry may be 
restored during inflation, and so, neither isocurvature fluctuation nor non-Gaussianity 
will arise. Also we show r = 10~^ by the green dash-dotted line. 



isocurvature perturbation is dominantly given by the linear Gaussian part S ^ S^^\ and 
hence the three point function is estimated as {SSS) ^ fs^s ^ /nl^^^- Thus we obtain 



On the other hand when 5a ^ FaO^ The non-Gaussian part dominates the isocurvature 
perturbation S ^ fsS^^^'^-^ giving a three point function as {SSS) ^ Pg^^ ^ InlP^- As a 
result we obtain 



In Figs, m [5] and [6l the non-linearlity parameter /nl is shown on H[^{-9 plane for 
Fa = 10^^ GeV, 10^^ GeV, and 10^^ GeV. It is seen that an observable amount of non- 
Gaussianity (say, /nl ^ 0(10)) is generated near the isocurvature constraint. However, 
notice that the relevant quantity from CMB observations is /nl^, not /nl, and the relation 
between them is given in Sec. [31 In particular it has been shown that can be 100 
times larger than /nl, depending on the observed scale. Thus it may be possible that 
isocurvature fluctuation is probed only through its non-Gaussian imprints on CMB. 




(71) 




(72) 
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Figure 5: Same as Fig. [4] with Fa = 10^^ GeV. The upper shaded region is excluded from 
the axion overproduction. 

5 Conclusions and Discussion 

In this paper we have investigated a possibihty that large non-Gaussianity is generated 
by isocurvature fluctuations. One interesting feature of this scenario is that the bispec- 
trum and the power spectrum of the CMB temperature fluctuations exhibit characteristic 
scale dependence. In particular, the eflFective non-linearity parameter f^^ is signiflcantly 
enhanced at large scales, compared to the adiabatic case. Furthermore, our results indi- 
cate that large non-Gaussianity may be accompanied with an observable fraction of the 
isocurvature perturbation. If future observations conflrm both large non-Gaussianity and 
a certain amount of isocurvature fluctuation component, our scenario will become very 
attractive. As a concrete example, we have shown that the axion can naturally induce 
such isocurvature perturbation in the CDM sector, leading to large non-Gaussianity. If 
the axion is indeed responsible for the large non-Gaussianity hinted by the current obser- 
vations, the inflationary scale should be in the range of 0(10^ — 10^^) GeV. This opens 
up an interesting possibility that the axion can be probed through its non-Gaussianity 
contribution to the CMB temperature fluctuation, even if the energy density of the axion 
today is negligible compared to the total dark matter abundance U- 

Although we have restricted ourselves to the CDM isocurvature perturbation in this 
paper, the baryonic isocurvature perturbation can also generate large non-Gaussianity in 

^ There may be an anthropic reason that forces the axion abundance to take such value [37] . 
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Figure 6: Same as Fig. H] with Fa = 10^^ GeV. The upper shaded region is excluded from 
the axion overproduction. 



a similar fashion [38] . There is indeed a scenario using the flat direction with the baryon 
number to generate baryonic isocurvature perturbations [391 1301 \^T\ . 

We have assumed that the CDM isocurvature perturbation comes from the fluctuations 
in the CDM sector. However, this may not be the case, and it can arise from the radiation. 
That is to say, the CDM isocurvature perturbation is given by S' = 3(Ccdm — Cr) = —^Cr 
in some cases, including the curvaton/ungaussiton scenario. In this case the sign of the 
non-linearity parameter fs becomes negative, which is expected to leave distinct features 
on the bispectrum of the temperature fluctuation. One may be able to use the features 
to probe into the origin of the isocurvature perturbation [12] . 

In general, the isocurvature perturbation mainly aflPects the large scale temperature 
anisotropy, while its effect on the density perturbation is weaker than the adiabatic one. 
Thus if the non-Gaussianity is truly sourced by the isocurvature perturbation, it can be 
seen only in the CMB observations, and we should have null detection of non-Gaussianity 
from the analyses using the matter power spectra. Hopefully, future cosmological obser- 
vations will enable us to establish or refute the existence of large non-Gaussianity. If it is 
there, we may be able to distinguish the origin, adiabatic or isocurvature. Undoubtedly, 
once detected, it will provide us with useful information on the early universe and the 
high energy physics. 
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